Radiative corrections and quantum gates in molecular systems 
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We propose a method for quantum information processing using molecules coupled to an external 
laser field. This utilizes molecular interactions, control of the external field and an effective energy 
shift of the doubly-excited state of two coupled molecules. Such a level shift has been seen in the two- 
photon resonance experiments recently reported in Ref. Here we show that this can be explained 
in terms of the QED Lamb shift. We quantify the performance of the proposed quantum logic gates 
in the presence of dissipative mechanisms. The unitary transformations required for performing one- 
and two-qubit operations can be implemented with present day molecular technology. The proposed 
techniques can also be applied to coupled quantum dot and biomolecular systems. 

PACS numbers: 03.67.Lx, 42.50.Fx, 32.80.Qk, 33.50.Dq 
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There is currently enormous interest in the field of 
quantum information processing (QIP). Potentially use- 
ful applications range from speeding up computations 
such as factoring, which would require tens of thousands 
of qubits, through searching and quantum simulations, 
which could be useful with a hundred or less, down to 
teleportation and quantum repetition of communication, 
which are effectively just few-qubit QIP. There have been 
some successful demonstrations of few qubit manipula- 
tions, and there exists a vast and still-growing range of 
proposals for realising QIP. The way forward is still very 
open — the key routes to few-qubit and large-scale QIP, 
which could well differ, have yet to be identified. Of par- 
ticular interest are new proposals that are based on al- 
ready observed phenomena, or that are clearly amenable 
to experimental investigation with current technology. 

In this Letter we show how to use small quantum net- 
works of coupled molecules to implement gate operations 
that enable universal QIP. A crucial ingredient to this is 
an effective energy shift of the doubly-excited state of two 
coupled molecules. We show that this can be explained 
in terms of the Lamb shifts of the levels, and we discuss 
how such a system can be tailored, with current technol- 
ogy, to realise a universal gate-set. We quantify the gate 
performance, and calculate molecular resonance fluores- 
cence spectra for comparison with experiment. Although 
the physical origin of the relevant energy level shift will 
be (potentially) different, it is clear that the approach 
can also be applied to other systems, such as coupled 
self-assembled quantum dots or biomolecules. 

We commence by describing our set-up. Consider 
two two-level closely spaced molecules fixed at positions 
Ti, which have excitation frequencies LUi, i = 1,2, and 
are separated by the vector r^. The dipole-coupled 
molecules are embedded in a dispersive medium of re- 
fraction index n and interact with the quantized radiation 
field and with an external coherent driving field of fre- 



quency ujl. Evidence of such dipole-dipole (d-d) coupling 
and the generation of sub and super-radiant states has 
recently been experimentally reported in a set-up that 
involves single terrylene molecules in a para-terphenyl 
crystal 0. We denote the ground and excited molecu- 
lar states by \gi) = \0i), and \ei) = with associated 
transition dipole moments d^, and corresponding spon- 
taneous emission rates P^. The strength of the molecule- 
laser coupling is given hy £i = — • E^, with E, being 
the amplitude of the coherent driving field at r^. The 
coupled molecules can be represented as a single four- 
level system in SU(4); we write the system interaction 
Hamiltonian in the computational basis of direct product 
states |z) (g) |j) ~ 0, 1) as 

= [i?^''+i?i'Vi^'^+i''^®[^o'^+^fV^i2 , (1) 

where hVi2 denotes the d-d interaction energy, and 
the free, laser, and interaction Hamiltonians are given 
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y ~ V) 01 denoting the transition operators, and 



by - 22 ~ ^ 11 

= ^^12 [Pi!^ ® PI2I + m2[Pl2 - ^ 21 

with Pi 

\i) and |j) being the eigenvectors spanning the Hilbert 
space. One can then map the transition operators into 
the pseudo-spin-1/2 dipole raising and lowering opera- 
tors S/, \ — -All , »s', — , — 



S^S^ = Pii , thus generating the corresponding opera- 
tor algebra in SU(2)®SU(2). We allow for a detuning 
Ai = LOi — lol between the incident field and the molec- 
ular transitions. We also incorporate a shift Ag of the 
doubly-excited state in our model — we shall show that 
this shift is an important feature of a robust molecular 
quantum gate. In our work we encode two qubits into 
this two-molecule system. Related work on such sys- 
tems, where a single qubit is encoded into a molecular 
dimer, is discussed in reference 0- In the following, we 
make use of the electric-dipole approximation fcori2 <C 1 
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(cfco = 27rt'o = 2tt{i/i + 1^2) /2) and the detuning parame- 
ters A_ = uji — UJ2, A+ = ll>i + UJ2 — '2ujl HI' 

The eigenstate coefficients of Hamiltonian are de- 
picted in Fig. n in terms of the different control param- 
eters of the model. We have used some of the experi- 
mental data reported in and some other input values 
that reveal interesting scenarios. The eigenenergy spec- 
tra (not shown) indicate that indeed the energy differ- 
ence between the entangled states shown in Fig.^a), for 
A_ = 2320 MHz (vertical line), is about 3 GHz, in agree- 
ment with the spectral results of However, notice that 
as the laser coupling is increased, there is a more signifi- 
cant contribution from all the basis states to the system 
eigenstates, thus producing (e.g., for A_ = 2320MHz) 
the superpositions aoo |00) -|-aoi |01) +aio |10) +aii |11) 
with weight coefficients as illustrated in Fig. ^a). 
This means that for strong laser pumping, aoo and an ^ 
in contrast to the entangled sub- and super-radiant 
states reported in 0. As shown in Fig.QJ such entangled 
states l^a) and \^s) ^ can indeed be created but for 
a different parameter window (e.g., for A_/Vi2 > 500). 
Figure n also shows that, for a fixed V12, the bigger the 
energy difference uji — UJ2, the easier one can achieve a 
'disentangling effect', where the system's eigenstates can 
become, to a very good approximation, equal to the com- 
putational basis states. This is illustrated as a function 
of the detuning parameter in Figs.^^c) and (d). 
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FIG. 1: Eigenstate coefRcients Uij of the coupled system in 
the absence and presence of laser pumping as a function of 
A_, and detuning A+/2. The weights of the basis states 
of the 4-level system have been plotted for V12 = 950 MHz 
(graphs (a) and (c)) and for V12 = 10 MHz ((b) and (d)). 
A_ = 2320MHz in graphs (c) and (d). Notice the difference 
in eigenstates behaviour for A_ ~ V12, and A-/V12 3> 1 0]. 

As a result of the dipole-dipole interaction between 
molecules, the system exhibits resonant cooperative ab- 
sorption of light 13, thus allowing simultaneous excita- 
tion by absorption of two-photons; we represent the spon- 
taneous decay of such excitations by the rate = r2i. 



It is not difficult to calculate the general expressions for 
the d-d interaction strength and the incoherent decay 
rate 0. In the near field (ri2 ^ \l) the inter- 
action energy and the incoherent decay rate are given 
by V12 = ^^^[di ■ d2 - 3(di • fi2)(d2 • fi2)], and 
ri2 — VI1I2 di-d2, and hence the maximum (minimum) 
strengths are obtained for parallel (perpendicular) dipole 
moments, z = nfcori2, F^ = = nu;f\\d\\^ / (Scohc^) , 

where Ai2-i (i = 1, 2) are the Einstein coefficients in vac- 
uum, and di and ri2 are unit vectors along and along 
ri2. The presence of the dielectric modifies the molecules 
radiative decay rates which in turn implies a shift in the 
molecular energy levels. The vacuum frequencies are now 
modified and hence the phase velocity of light waves in 
the dielectric medium is given by c/n{uj) = uj/k{uj). 

We describe the full molecular quantum dynamics by 
solving the master equation for the Lindblad operator [3] , 
thus taking into account all the possible excitation de- 
cay channels F^, and Fy. In particular, Ref. Q re- 
ports an energy shift in the fluorescence spectrum as- 
sociated with the upper level of A-|_/2 ~ —160 MHz 
(see the central resonance of Fig.|2Ia)), where the steady- 
state density matrix element Pee.ee = (eie2|p|eie2) dis- 
plays the observed shift (A_|-/2 = —160 MHz, verti- 
cal dashed line). This effective shift can be attributed 
to the Lamb shifts of the three states \^a), |^s) and 
Y^c)- Decay channel interference modifies the decay 
rates of the first two states [J (and hence their Lamb 
shifts) away from the uncoupled values. As a conse- 
quence, the doubly-excited state is effectively displaced 
from where it would be expected, given the energies of 
the lower two levels. The Lamb shifts can be calculated 
by taking into account self- interactions and by renor- 
malizing to remove quadratic divergences |^ |3| from 

A. = -iS#/"H[(^+^o-; + (^-c^»r']rfc.. A 

detailed calculation for the experimental set-up reported 
in shows that A^ = A^ - A^ - Aa ~ -178 MHz, a 
result that is in good agreement with the observed value, 
as is shown by the inset of Fig. |21 where the fluorescence 
around this excitation energy is plotted for both the the- 
oretical (solid line) and the measured (dashed line) values 
as a function of the laser detuning A_|-/2. As we show 
below, the A,; shift is clearly useful for QIP; even more 
so because it can be tailored by means of an appropri- 
ate engineering of the structural and optical properties of 
the molecular set-up. Different strengths A^ are shown 
in Fig. |2Jb) where the expected fluorescence spectra is 
plotted as a function of the laser detuning A_|_/2. The 
result that additional energy shifts appear in the spectra 
is viewed here as a meaningful resource for performing 
molecular two-qubit logic gates and therefore for quan- 
tum computing purposes. Although we have focused on 
the experimental data reported in 0], other parameter 
windows and different molecular systems can certainly 
be explored in the light of our results. 
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The splittings in the energy spectrum of the inter- 
acting system can be used to induce a conditional dy- 
namics between the two subsystems (e.g., a laser with 
Al ~ 578 nm and linewidth ~ 1 MHz can be suit- 
able 0). Figure H] displays the eigenstate coefficients in 
the absence and presence of laser pumping for different 
interaction strengths V\2- For li = 0, the energy sep- 
aration between 1*^) and |*a) is \J ^- + 4] and 
hence molecules of equal transition energies (A_ = 0) 
have Ci = C2 = 1 and the corresponding symmetric and 
antisymmetric states become the maximally entangled 
states (MES) -^{\e\g-2) ± \g\e-i))-, as can also be seen 
in Figs, da) and (b). Otherwise (wi ^ ^2), the degree 
of molecular entanglement can be manipulated by con- 
trolling the ratio Vi2/A_. This can be experimentally 
done by means of applying a differential Stark shift via 
an external inhomogeneous electric field to select differ- 
ent values of the detuning A_ or by directly tuning the 
dipole-dipole interaction by means of changing the sepa- 
ration ri2 and/or by using a medium with different dis- 
persive properties. This means that a rich spectra of en- 
tangled symmetric and antisymmetric states can be gen- 
erated by modifying the eigenstate coefficients. It follows 
from that ai = ^(F + l)/2r, = ^/{Y -1)/2Y 
{Y = X/A_), then, if the system's ratio Vi2/A_ < 1, 
a2 « Fi2/A_ {ci/^/2 = ai): the bigger the difference 
in the transition energies (compared to V12) the more 
disentangled the intermediate states become. 

The effective shift A^ described above implies that 
the resonant frequency for transitions between the ba- 
sis states \gi) and jci) of one qubit depends on the state 
of the neighbouring qubit: the energy difference between 
states I'^a) and that of is different from the one be- 
tween states and \"^s) due to Aj and V12 0. Thus, in 
the scenario Vi2/A^ ^ 1, we can naturally construct the 
two-qubit CNOT gate J7cnot: \m) ^ \n) 1-^ |to) (g) |m ffi n). 
Hence, the logic operation |ei) \g2) i-+ |ei) [62) can be 
achieved by illuminating the input qubit system with a 
TT-pulse of energy SI12 = lu2 ~ S + (say, a Trnij-pulse), 
where 6 = V12/A.-. This is illustrated, for experimen- 
tally attainable values 0, in Fig.|3Ia) where the system's 
initial state \e1g2) is induced to evolve onto the final state 
\e1e2) via the application of a 7rnj2-pulse, thus realizing a 
quantum CNOT gate in 1.25 ns. In the absence of dissipa- 
tive channels the gate operation exhibits a perfect fidelity 
(see below) , whereas the inclusion of realistic decay rates 
shows the appearance of small weight contributions aris- 
ing from the populations Pgg,gg and Pge,ge (curves at the 
bottom of Fig.|3fa)). This, however, can be further con- 
trolled by means of choosing different Aj, ¥12/ A^, and 
laser strengths to improve the gate fidelity. Conversely, if 
the role of the control qubit is to be performed by the sec- 
ond qubit, the gate operation \gi) |e2) |ei) |e2) can be 
realized via a Trjiai -pulse of frequency VL21 — uJi-^ 5 + A.^. 
The energy selectivity of such a CNOT gate is thus deter- 



mined by A_, Ag and V12, and typical operation time 
scales can be gathered from Figs. |5] and for differ- 
ent parameter windows. The complimentary unitaries 
in SU(2) required to build a universal set of gates can 
be achieved by inducing appropriate Rabi oscillations 
of the corresponding dipole operator expectation value 
(d) = (Pi2di2 -I- f2id2i). Thus, single qubit control can 
be performed by tailoring the energy and length of ap- 
plied laser pulses to induce distinct Rabi oscillations. 
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FIG. 2: Steady-state density matrix elements as a function of 
the laser detuning parameter A+/2. The occupation proba- 
bilities give the structure of the resonance fluorescence spec- 
trum, as is shown in a) for different excitation intensities li, 
and in b) for different effective shifts A^; i!i = il, Ti = 50 MHz, 
Fia = 9MHz, A_/2 = 1160 MHz, and V12 = 950 MHz. 



The scenario considered previously exploits the diag- 
onal shifts Aj by means of focusing on a particular pa- 
rameter window. However, even for arbitrary V12/A.- 
ratios one can still perform entangling gates that enable 
universal QC. By writing the interaction Hamiltonian of 
Eq. (P) as H12 = hVi2[crx <^ cr^ + cTy ® ay] = H^y, one 
can define an optimal entangling gate Uxy for the d-d 
interaction by simply letting the system evolve under 
Hxy for a time txy = tt/AVi2. By means of the sclf- 
adjoint operator Hp J2if^i^i = exp(-ii7p:) 

provides, up to LU, a decomposition for any two-qubit 
gate and thus the d-d interaction has associated Uxy'- 
\m) ®Jn) I— > \n) (g) \m) as its 'natural' entangling 

gate Isj. This gate, equivalent to a combination of CNOT 
and SWAP, plus LU, are a universal gate-set. This means 
that, provided one has access to performing sufficiently 
fast LU operations, one could always exploit the d-d in- 
teraction for universal QC. Typical switching times for 
the experimental data analyzed here can be in the sub- 
picosecond time scale (~ 1/1000 the fastest decay rate). 

We calculate the steady-state density matrix elements 
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in the computational basis as a function of the laser de- 
tuning A+/2. Fig. |5fa) shows a strong enhancement of 
the two-photon resonance (TPR) as the strength of the 
laser coupling is increased (for V12 = 950 MHz). As the 
d-d coupling determines the size of the TPR, the weaker 
the interaction the more difficult the detection of the 
TPR fluorescence peak. Fig. ^h) shows effectively that 
an increase in has associated with it a further TPR 
shift that could then be observed in the fluorescence spec- 
trum, as predicted above. The asymmetric shape of the 
resonance lines arises from the destructive and construc- 
tive interference between the system's excitation decay 
channels ^ \^a) ^ l^^e), and 1^-^) |*,), 

which in turn is also reflected in the observed linewidth 
difference between the left and right peaks (Fig. EJ, a 
hallmark of Dicke's subradiance and superradiance. 
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FIG. 3: (a) CNOT gate operation: a yrnj^j -pulse is applied in 
order to generate the conditional logic \e\g2) ^ |eie2). The 
corresponding populations are shown as a function of time in 
the absence and presence of dissipative channels. A+/2 = 
-1319 MHz, A_ = -2320 MHz, h = = 200 MHz, A, = 
160 MHz, and V12 = 50 MHz. (b) Fidelity (see text) in the 
generation of the MES -^(|5i52)-l-|eie2)). The bar signals the 
end of the applied ttq^j -pulse and indicates a measure of the 
actual entangled state achieved in the absence and presence 
of dissipation. Input parameters are as before except for = 
50 MHz. Decay rates are given in 10^ s~^ units. 



The quality of the gates discussed here depends on 
the ability to control the relevant structural and exter- 
nal parameters. An indication of gate performance is the 
quality of the maximally entangled states they produce. 
Due to decoherence, any attempt to generate a MES 
IV'm) in practice creates a mixed state p. In Fig. Ofb) 
we plot the corresponding fidelity T = ^ (V'm|p|'0m) for 
\'il)m) = :^(l5i52) + |eie2)). The initial state is taken to 
be -^(|.9i) + |ei)) 152) (this can be achieved by starting in 
the system's ground state and then applying a 7r/2 pulse 
at energy wi), and hence IfAm) can be generated by apply- 
ing a 7i"ni2-pulse. Figure |3fb) gives the fidelity between 
IV'm) and p in the presence of dissipative interactions. 
Clearly, zero decay gives a perfect fidelity [T = 1), but 
the action of dissipative processes (as described in 0) 
gives T 1, as shown in Fig. |3Ib) by the bar signaling 



the end of the applied tthj^ -pulse. We have considered the 
rather non-optimistic ratio ~ 2 x 10~^ but smaller ra- 
tios could indeed be experimentally achievable and hence 
improved fidelities can be obtained in conjunction with 
tailored A^, and laser field values. Our numerical ex- 
periments consider a setup and energy scales that have 
been measured in the laboratory Q thus providing a sig- 
nificant step towards proof-of-principle experiments that 
can be undertaken with currently available technology. 

Although of a different physical origin, we stress that 
the main ideas discussed here can also be applied in the 
context of the technologically irnportant self-assembled 
coupled quantum dot molecules [liJ, and biomolecu- 
lar systems (l2| . Summarizing, we have given a detailed 
prescription for molecular conditional quantum dynam- 
ics that can be implemented with modest experimental 
requirements by tailoring control parameters such as ex- 
citation frequencies, spontaneous decay rates, interaction 
strengths, and external laser pumping. Furthermore, we 
have explained the quantum mechanical origin of the 
TPR energy-level shift and we have highlighted its rele- 
vance to QIP protocols. We have shown that the perfor- 
mance of a molecular quantum gate can be adjusted by 
varying the value of this shift. 
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